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Abstract 

A translational gauge approach of the Einstein type is proposed for obtaining the 
stresses that are due to non-singular screw dislocation. The stress distribution of sec- 
ond order around the screw dislocation is classically known for the hollow circular 
cylinder with traction- free external and internal boundaries. The inner boundary sur- 
rounds the dislocation's core, which is not captured by the conventional solution. The 
present gauge approach enables us to continue the classically known quadratic stresses 
inside the core. The gauge equation is chosen in the Hilbert-Einstein form, and it 
plays the role of non-conventional incompatibility law. The stress function method 
is used, and it leads to the modified stress potential given by two constituents: the 
conventional one, say, the 'background' and a short-ranged gauge contribution. The 
latter just causes additional stresses, which are localized. The asymptotic properties of 
the resulting stresses are studied. Since the gauge contributions are short-ranged, the 
background stress field dominates sufficiently far from the core. The outer cylinder's 
boundary is traction- free. At sufficiently moderate distances, the second order stresses 
acquire regular continuation within the core region, and the cut-off at the core does 
not occur. Expressions for the asymptotically far stresses provide self-consistently new 
length scales dependent on the elastic parameters. These lengths could characterize 
an exteriority of the dislocation core region. 



PACS: 03.50.-z, 11.15.-q, 61.72.Bb, 61.72.Lk 



1 Introduction 



A model of non-singular screw dislocation is proposed in the present paper. The description 
is based on the translational gauge approach of the Einstein type. The aim of the paper is 
to investigate the corresponding stress problem in the second order approximation. 

Solutions to the stress problem for the edge and screw dislocations by means of the 
theoretical elasticity arc known both in the linear and quadratic approximations (see [1]- 
[4]). The stress distributions of the first order are singular on the dislocation lines. Beyond 
the scope of the linear elasticity, the dislocations have been conventionally considered in [5]- 
[8] (the non- linear approaches are reviewed in [9] and [10]). The stress solutions provided 
by [5]- [8] are valid for the hollow cylinders subject to the traction-free conditions on the 
external and internal boundaries. The corresponding internal radii, as the dislocation core 
boundaries, remain unspecified. Note that 'core' is rather a crystalline notion since the 
theoretical elasticity does not have an appropriate length scale. 

The relevance of the second order effects in the theoretical elasticity to the lattice mod- 
eUing of the dislocations has been discussed in [9]. It is not sufficient to use only the 
linear elasticity for the description of strongly distorted regions of the dislocation cores (see 
also [3], [4], [9], and [10]). Self-consistency of appropriate boundary- value problems at the 
dislocation cores also seems theoretically challenging. 

The Lagrangian field theory based on the gauge group r(3)E) S0{3) has been proposed 
in [11] and [12] as a non-conventional approach to the defects in the elastic continua. The 
Maxwell-type Lagrangian has been proposed to govern the translational gauge potentials. 
The ordinary screw dislocation, as an asymptotic configuration, is allowed for in [11] and [12]. 
Owing to the relevance of higher approximations, it has been further attempted in [13] 
to use [11] for obtaining the corresponding stresses of the second order around the screw 
dislocation. 

An important modification has later been proposed in [14] for the translational gauge 
equation advanced in [11] and [12]. The modification is concerned with the corresponding 
driving source. In its new form, it is given by the difference between the stress tensor of the 
model and the stress due to a classical background dislocation. In this particular case, the 
choice of the background corresponds to the screw dislocation. The linear approach of [14] 
results in a modified screw dislocation, which is equivalent, for sufficiently large distances, 
to the ordinary one. It is crucial that additional short-ranged stresses, which remove the 
classical stress singularities, arise. However, it is problematic [15] to obtain analogously a 
non-singular edge dislocation by means of the translational Lagrangian [11] alone. 

Further, the translational gauge approach to the elimination of the dislocation singulari- 
ties has been developed in [15]- [19]. The driving source of the gauge equations in [15]- [19] 
is the same as in [14], i.e., it is given in the form of the difference of two stresses. The same 
modified screw dislocation as in [14] is allowed for in [15] and [17] as well. 

Non-local elasticity should be also recalled as a non-gauge opportunity for extending the 
conventional elasticity framework. Similarly to [14], the singularity-less screw dislocation 
has already been independently obtained in [20] (see also [21]) by means of the non-local 
elasticity. Moreover, the first strain gradient elasticity [22]- [27], as well as the generalized 
elasticity [28]- [30], also admits non-singular solutions for the dislocations. 

Let us turn again to the gauge approach [13], where the method of stress functions [5], [7] 
is used. In principle, the corresponding Kroner ansatz for the screw dislocation, in the second 
order, is the same as for the edge dislocation, in the first order. However, according to [15], 
the use of the Maxwell-type translational Lagrangian does not allow for the Airy stress 
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function with a correct numerical coefficient. Although this Lagrangian is successful for the 
screw dislocation in the ffist order, its usage beyond the scope of the linear approximation 
looks doubtful. Moreover, the elastic energy in [13] does not contain the terms of third order, 
which are accounted for in [5], [6], and [7]. There is no continuation of the stresses within 
the core region. Therefore, the experience of [13] looks, to a certain extent, incomplete. 

In principle, the translational gauge approaches discussed in refs. [15]- [19] look promis- 
ing for study in the second order. Indeed, the length scales are "generated" in these models, 
and the singularities are smoothed out within compact core regions. The active interest in 
modelling the defects and their cores, as well as the importance of the higher approxima- 
tions, stimulate us to investigate the gauge models further in order to gain more realistic 
descriptions of the core regions. 

Specifically, ref. [15] is based on the Hilbert-Einstein gauge Lagrangian. Two solutions 
are found therein, which modify, for short distances, the Airy and Prandtl dislocational stress 
functions. The stress singularities of both classical dislocations are thus removed. However, 
certain components of the modified stresses in the "edge" case do not behave properly at 
infinity. 

The approach used in [16]- [19] is based on the Lagrangian quadratic in the r(3)-gauge 
field strength. This quadratic form is based on three gauge-material constant^ and therefore 
is different from that in [14]. The inclusion of the terms corresponding to the energy of the 
rotation gradients enables one to obtain the modified edge dislocation with an appropriate 
large distance behaviour (see [19]). The model from [19] contains, in general, ten material 
constants. 

The model [15] looks attractive for study in the quadratic approximation. Indeed, the 
Kroner ansatz for it can be specified, in each order, appropriately. The gauge approach, 
say, [16]- [19], would require a rather sophisticated choice of the gauge parameters for dif- 
ferent orders. The Hilbert-Einsteinian approach implies a single gauge parameter and looks 
preferable because of the essential non-linearity of descriptions based on it. Besides, the 
corresponding Einstein tensor is important in the conventional dislocation theory as well. 

In the present paper, the model [15] is investigated along the line of [5], [6], and [7]. The 
background defect is given by the straight screw dislocation in the infinite circular cylinder. 
Its stress field is taken up to the second-order terms considered in [7]. The Einstein-type 
gauge equation plays the role of unconventional incompatibility law. The stress function 
approach is used. The stress function of second order is obtained. It is equal to the sum of the 
background and non-conventional parts. The latter results in a continuation of the stresses 
within the core region, and the cut-off at the core's boundary is avoided. The asymptotic 
properties of the second order stresses are studied. New length scales of the "gauge" origin 
are provided, which could characterize an exterior structure of the dislocation core. 

We do not discuss in detail a complicated theme such as the dislocation cores (see [3], 
[4], [9], [10], [32], and [33]). For instance, [34] and [35] are the first attempts to incorpo- 
rate the lattice periodicity and to display the finiteness of the cores. Early attempts at 
treating the screw dislocations and their cores atomistically are given in [36], [37], and [38]. 
A singularity-less screw dislocation has also been obtained in [39] by means of the quasi- 
continuum approach (see [40] as well). In their turn, the gauge fields are important in modern 
condensed matter physics (see [41], [42], and [43]). Certain applications of [11] and [12] can 

^For a special choice of the parameters, the gauge Lagrangian in [16]~ [19] is equivalent to the Lagrangians 
either in [14] or in [15]. Note that the most general eight-parameter three-dimensional geometric Lagrangian 
is known [31] , which includes, in addition to the Hilbert-Einstein term, the terms quadratic in the components 
of the differential-geometric torsion and curvature. 
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be found in [44]- [49]. 

The paper consists of six sections. Section 1 is introductory (see also [14]- [19]). Section 
2 outlines the Einstein-type gauge equation. Section 3 specifics the pcrturbativc approach. 
The gauge equation is solved, and the modified stress potential is obtained in Section 4. 
The stress fields and their asymptotics are investigated in Section 5. Section 6 closes the 
paper. Details of the calculation are given in Appendices A and B. The present paper is a 
refined and improved (concerning the large distance behaviour of the stresses) version of [50] . 
Bold-faced letters are reserved for second-rank tensors, and appropriate indices can easily 
be restored. 



2 The Einstein-type gauge equation 

The Einsteinian gauge equation is outlined in the present section. For more on the appro- 
priate differential geometry, one should refer to [10], [42], and [51]. The Eulerian picture 
is chosen in the present approach instead of the Lagrangian one accepted in [15]. The ini- 
tial and final states of the dislocated body are referred to the coordinate systems {x*} and 
{x"}, respectively. The corresponding squared length elements are expressed as g.ijdx'^dx^ 
and Tjabdx^'dx'^ (throughout the paper, the indices repeated imply summation) . Let us define 
the frame components e" by the relation di — e"5a (henceforth the partial derivatives d/dx'^ 
arc denoted by di). The co-frame components are given by the one-form dx* = S^dx"-. 
The components and their duals e" are orthogonal: e"£^' = 5^, e'^^Sj — Sj. For the 
metric r]ab we get r]ab = gij^a^b ■ 

The Eulerian strain tensor is referred to the deformed state, and it measures a deviation 
between the final and initial configurations [10]. Let us map an initial state to the deformed 
one: x i — > ^(x). Then the difference 

VabdCdi' - gijdx'dx' = 2eabdCdi\ (2.1) 

where 

2e„b = r]ab - gab: gab = gij^^B^ , (2.2) 

defines the Eulerian strain tensor {e)ab = Gab- The metric Qab in (2.2) is called the Cauchy 
deformation tensor. Here are given by the one-form dx* = B^d^"' as follows [10]: 

^ ^ = - Va , (2.3) 

provided that the displacement ^ — x is expanded as — x'^ = u°'S^. The covariant derivative 
Va in (2.3) is defined by the requirement that the components are covariantly constant: 

Va£b'^da£b'-[^^] ^; = 0, (2.4) 

where | | is the Christoffel symbol of second kind. The metric rjab — ^a^bi is also 

covariantly constant because of (2.4). In turn, (2.4) allows one to express the Christoffel 
symbols in terms of the metric rjab (see [42]). 
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To implement the gauging of the group T(3), we extend the co-frame components BJ^ 
(2.3) by means of the so-called compensating fields: 



'a 



1^ - = ^: - {^c va u'' + . (2.5) 



The entries ^p^ are the translational gauge potentials, which behave under the local trans- 
formations — > -|- r]'^{x) as follows: 

9x* dx^ { drf \ 



Qta ' Qta \ 3 ^ Q^j ) ' 

(2.6) 

fa ' ^ d^dx^' 

The shifts (2.6) ensure the gauge invariance of i3J (2.5) (and of the definitions (2.2), as well). 
Motivation of (2.5) should be traced back to T(3) E) 5'0(3)-gauging and to the corresponding 

Cartan structure equations [11] (see also [52] and [53]). The occurrence of the non-trivial 
core region for the screw dislocation in question (i.e., the appearance of the short-ranged 
stresses) should be associated with the gauge variables ip^. The Eulerian strain (2.2) takes 
the gauged form: 

in) (ri) Av) \ fiv) \ 

'^eab ^VaUb + fab + V b Ua + fba " [V a Ub + fab J [V b Ua + ^ba) ■ (2-7) 

When cp is zero, (2.7) is reduced to the conventionally-looking (background) strain tensor in 
the Eulerian picture [10]. 

The differential-geometric approach of the present paper is the same as in [15], i.e., the 
so-called teleparallel framework (see [52] for a quick remind) is adopted as the most suitable 
for describing the dislocations [51]. Let us consider the Riemann-Christoffel curvature tensor 

Rd. 
abc ■ 

where the Christoffel symbols are calculated for the metric Qab (2.2) (the subscript 'g') given 

. , , (g) (g) 

by B^ (2.5). The metric gab is covariantly constant, i.e., Va Qbc = is fulfilled, where Va 

iv) 

is defined similar to (2.4). We obtain from Va Vbc — the following relation between the 
Christoffel symbols: 

c 1 f c 



obi lab' -2-=- ■ P.9.1) 

/(»?) iv) (»?) \ ^ ^ 

2eab'' = g [Va ebe+ Vb eae- Ve Babj ■ (2.9.2) 

where (2.2) is taken into account [10]. 

Let us define the Einstein tensor G^^ = - £^°'^£^'^Rabcd, where S^^'^ is the Levi-Civita 

tensor defined by means of the metric rjab- Then, the Einstein-type gauge equation [15] takes 
the form 



Here £ is a constant factor at the Hilbert-Einstein gauge Lagrangian. The geometry of the 

(v) , 

deformed state is Euchdean in the Eulerian approach. The corresponding curvature Rabc 

is zero for the metric rjab- We substitute (2.9.1) into (2.8) and use the vanishing of Kabc^- 
Then, (2.10) is re-expressed : 

(r;) (ri) 1 

_^eab^fcd ^ _ (^e/ _ (^^je/) ^ 2E^-'£^'^''eadl e^c' • (2.11) 

The differential-geometric torsion, as an independent degree of freedom, is not considered. 
It is just the Riemann-Christoffel curvature that is subject to the gauge equation (2.10). 
Therefore, (2.11) governs only the variables related to the metric, i.e., the corresponding 
strains. 

The variational derivation of (2.10) can be discussed along the hne of [15]. The right- 
hand side of (2.10) is given by (2£)~^(cr — cTbg), where (Tbg implies the background stress 
tensor. The deviation of the stress tensor of the model from cr bg plays the role of the driving 
source in the gauge equation. The parameter ^ characterizes an energy scale intrinsic to the 
gauge field ip: it makes the driving source dimensionless. In the present paper, cr^g implies 
the stress field of the straight screw dislocation lying along an infinite cyhndric body. Both 
cr and (Tbg respect the equilibrium equations: 

Vaf7"^ = 0, VaK)"' = 0. (2.12) 

A more detailed consideration of the gauge geometry behind the model in question should 
be done elsewhere, but several refs. should be mentioned in addition to those listed in [15]. 
For instance, useful indications concerning the translational gauge geometry can be found 
in [54] and [55]. A topological picture is proposed in [56], which includes the dislocations 
and extra-matter by means of the torsion and non-metricity, respectively. Moreover, [57] 
provides a further development of the geometric approach [31]. 



3 Specification of the gauge equation 
3.1 The stress function method 

We shall investigate the gauge equation (2.11) using the method of stress functions proposed 
in [5] and [7] for solving the internal stress problems in the incompatible elasticity. An 
exposition of [7] can be found in ref. [10] devoted to a review of the dislocation problems 
in non-linear elasticity. Certain details omitted below can be restored with the help of [7] 

and [10]. 

We shall use the stress function approach in a successive approximation form. Let us 
represent the strain and the stress tensors perturbatively: 

(1) , (2) (1) , (2) 

ef^ie-Fe, cr ^ cr + cr , (3.1) 

(2) (2) (1) (1) 

where e and cr are of second order smallness in comparison with e and cr. The repre- 
sentation (2.7) implies that each contribution in (3.1) consists of two parts: owing to the 
background and to the non-conventional origin. Following [7], we shall investigate only 
the stress problem. Finding the relationship between the non-conventional stresses and the 
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gauge potentials is beyond the scope of the present paper. Substituting (3.1) into (2.11) and 
(2.12), we obtain equations of the first {i — 1) and second {i — 2) orders: 

V-S^O, (3.2.1) 

H) 1 (i) 

Inc e = — ^ ir + Q . (3.2.2) 

The tensor notations [10] are used in (3.2). For instance, the left-hand side of (3.2.1) takes 
the form of the divergence of the tensor of second rank. The notation in (3.2.2) is: 



^Inc 



(^g)%i;'--(g,,)"', (3.3.2) 

J Uf (1) (1) , 

in) 

where 'Inc ' is the double curl incompatibility operator, V is the covariant derivative Va 
(= Va), and the indices are raised and lowered by means of the metric rj. 

Equations (3.2) look similar to the conventional equilibrium and incompatibility laws: see 
Eqs. (622) and (623) in [10]. The equilibrium equations of the first and second orders are 
given by (3.2.1). The gauge equations (3.2.2) play the role of non-conventional incompatibil- 

ity conditions. However, there exists a distinction that is due to {2i)~'^d a in (3.2.2). This 
term is just responsible for the short-ranged behaviour of the gauge parts of the resulting 

(1) 

stress functions. Moreover, Q in (3.3.3) is free from the torsion tensor in our approach. 
The elastic energy of an isotropic body is chosen in the Eulerian representation as follows: 

W{e) = jll{e) + khie) + l'll{e) + m'h{e) h{e) + n'h{e) , (3.4) 

where j = + A/2 and k — — 2/x are the elastic moduli of second order (A and ^ are the 
Lame constants), and V , m! , and n' are the elastic moduli of third order. Once W{e) has 

been chosen in the form (3.4), the constitutive law relates e with <t (i = 1, 2) as follows [7]: 

e = Cih((T)r] + C4 <r + * , (3.5) 

(0) 

where and 

*^ + C^I,(S))ri + C5/i(S) S +C,h{^S) . (3.6) 

Here, the /i,2,3(-) are the tensor invariants [58]. The numerical coefficients Ci and C4 are: 

Ci = -(2/.)"^^, C4=(2/.)-\ (3.7) 

where u = A/(2(A + fi)) is the Poisson ratio. The coefficients C2, C3, C5, and C7 can be 
expressed (see [7] and [10]) in terms of the elastic moduli of second and third order (3.4), 
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although this is not used below. More details on the relation of A, /i and m', n' with the 
elastic moduli of crystals can be found in [4] and [10]. 

The equilibrium equations (3.2.1) should be fulfilled with the help of the Kroner ansatz: 

<T = Inc X . (3.8) 

Substituting the constitutive relations (3.5) into (3.2.2) and using (3.8), we obtain equations 
for the components of the stress potential 

A A Xab +a(VaV;. - r^„;,A) A7i (g) + ((1 - a)V„V5 + a ?7„6 a) V^V^ X^d- 

(3.9) 

where 

(1) (2) (1) / 

S ab = 0, S (ab) = Q (ab) - InC * 

\ / (ab) 

(j) 

A is the Laplacian, and the S a are expressed by means of (3.3.2) and (3.8). In addition, 

(1) (1) 

Q and 'I' are given by (3.3.3) and (3.6), accordingly. The curly brackets around the indices 
imply symmetrization, and we use the dimensionless parameters i<? = and a (see [11] 
and [12]): 



3A + 2// 1 + 3A + 2// 1 + 1/ 

3.2 The gauge equations in the first and second orders 

Let us adjust (3.9) to the special problem in question. We replace the derivatives Va by the 

partial derivatives (9a, where are the coordinates in the final state, and assume 93 = 
(see [7] and [10]). Conventional notation is adopted for the components of the stress potential 
which are non-trivial [51]: 

{i) (i) (i) 

= ^2 Xi3 - <9i X23 , « = 1, 2 , 

(3.11) 



/ = X33, p= -dfi X22 Xii +23^2 X 
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The other components of x zero. The background stress tensor cTbg is also given by 
(3.8), but the corresponding stress potential is labelled appropriately: Xhg- 

3.2.1 The first order 

(1) 

In the first order, only the stress potential (3.11) is non-zero for the screw dislocation. 
The corresponding governing equation appears from (3.9) as follows [15]: 

(1) (1) 

{A-k'){<I> - (/.J =b5^'\x), (3.12) 
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(1) 



where 0bg= (~^/27r)logp is the background stress potential, and b is the Burgers vector 
length. The Burgers vector is parallel to the line of the screw dislocation, and its length is 
the perturbative expansion parameter. The axial symmetry of the final state stimulates the 
usage of the cylindric coordinates p, (p, and z instead of {x"} (p and if are chosen in the 



1 2 
X , X 



-plane and z 



Xs). We obtain from (3.12) the stress potential of the first order: 

(1) (1) 



0bg -fs, fs = {b/2TT)Ko{Kp) 



(3.13) 



where Kq is the modified Bessel function [59]. 

The stress a^z is non-trivial for the modified screw dislocation in the first order, and we 
express it, using (3.8) and (3.11): 



- pdp 



(1) 



hp 
2^ 



P 



'(l - KpKi{Kp)) 



(3.14) 



Equation (3.14) demonstrates a core region at p< k ^ : the gauge correction to the classical 
long-ranged law 1/p is exponentially small outside this region. Inside it, the law 1/p is 
replaced by another non-singular one. More detailed information on the numerical behaviour 
of (3.14) (including a treatment of in terms of interatomic spacing) can be found in [14] 
and [17]. In the present paper, it is assumed that b = 0{k~^). Solution (3.14) is in agreement 
with [14], [17] (the translational gauging) , [20] (the non-local elasticity) and [22], [26], and [27] 
(the strain gradient elasticity). This is because the Helmholtz-type governing equations 
similar to (3.12) are essential in the approaches mentioned. 

Certain analogies in the structure of the tensor laws governing the electrostatic of di- 
electrics and the elastostatic problems have been discussed in [60]. It looks hopeful that 
the short-ranged constituent of the first order solution presented here is also falling into the 
class of problems considered in [60]. In addition, within an independent framework though, 
strongly localized stress potentials turn out to be also responsible for the Debye-like screening 
effects in the dislocation arrangements [61]. 

3.2.2 The second order 

In the second order, from (3.9) we obtain the following gauge equations d: 





(1 - a)p + 


aAf - 


n'if- 


/bg 




(1 - a)p + 


aAf - 




/bg 




(1 - a)p + 


aAf - 




/bg 


1 - a) AAf + aAp - 




-Phg) 





2p 



2p 



Qii + 922^33 



Q22 + dl,^3s 



(3.15.1) 



2p 



Q33 - 5?2(^12 + ^21) + ^11^22 + dl^^u 



(3.15.2) 



Equations (3.15) describe the stress potentials of second order / and p (3.11). In addition, 

(2) (2) (2) 

there are equations to determine 0. But since (f)^„ is zero [7], we put consistently 0=0. 



^Henceforth, Q and * are used without the superscript. 
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Further, it is necessary to express Q in (3.3.3) and * in (3.6) in terms of the first order 
solution (3.13). Taking ^3 = into account and using tah^ (2.9.2), we obtain the following 
non-zero components: 

(A 6) 

Qii = Q22 = ^ , (3.16.1) 

/ , (i)x2 , (1) , (1) (A 0V 

In turn, the corresponding components of ^ are expressed standardly (sec [7] and [10]). 
Using (3.16.2), we obtain the combination that we are interested in: 



where 



= A^33 + (1 - 2At'C7) $ + 



, (3-17) 



/ (1)\2 (1) (1) 



(3.18) 



(1) 

Equations (3.16) and (3.17) look unconventionally owing to the presence of (A 0)^/4. 
It is just the dependence of Q on the torsion that results in the absence of the contribution 

(1) 

(A 0) /4 in the components Qii, Q22, and Qss obtained in [7]. Recall that the conventional 
incompatibility law requires the vanishing of the Einstein tensor calculated by means of the 
Riemann-Cartan curvature [51]. This latter includes (see [42]), in addition to the Ricmann- 
Christoffel part (2.8), a contribution owing to the torsion tensor. The latter is identified as 
the dislocation density. That is, the dislocation density contributes into the driving source 
of the incompatibility law. On the contrary, according to [14], [15], and [17], a stress field 
owing to the background defect is chosen first. Since the torsion is not considered as an 
independent variable, Qu, Q22, and Q33 acquire the form (3.16). 

Equations (3.15) allow for a correct limit to their classical version, which is respected by 
the background stress potentials of second order /bg and pbg- Let us pass to the classically- 
looking Qn, Q22, and Qu- 

Qn^dl^Q', Q22 = dlQ', Qu = -dl^Q' ■ (3.19) 

(2) 

A specific value of the constant Q arises from the requirement that a 33 averaged over the 
cylinder's cross-section is zero (sec [7] and [10]). Then, the equations that govern /bg and 
Pbg take in our notation the form (ac 0) 



AA/bg = k [A(^^f -I- Q') + (1 - a) (1 - 2/.2C,) 
(1 - a) pbg + a A/bg = -2/x (*^f + Q') , 



(3.20) 



where 

2/i 1 + ly , . 

k = — — = 2a^—. (3.21) 

l-2a ^ ^ 



10 



b '■^^ 

Moreover, and $bg are given by (3.18), provided that is replaced by cf)^^. 

Let us turn to (3.16). In the present approach, (A 0j is the square of the density 
distribution of the modified screw dislocation (see [14], [15], and [17]): 

(A^^y = (^±.^Ko{Kp)y. (3.22) 

For large distances, Kq{kp) decays exponentially, and (A j given by (3.22) can approxi- 
mately be replaced by zero. Therefore, the components (3.16) take, also approximately, their 

(1) 

conventional form. Then, three equations (3.15.1) reduce to a single one. However, A is 

not negligible for sufficiently moderate Kp, and some care inside the core region p < is 
required. In other words, this might be an indication for extension of the geometric frame- 
work by means of the differential-geometric torsion. But this would, in turn, imply that the 
teleparallel description is abandoned. 

Instead, we shall estabhsh an "effective" picture, which should be viewed as still contained 
within the teleparallel framework. However, this description is expected to incorporate 
certain features of the approach extended differential-geometrically. Indeed, an independent 
non-trivial contribution owing to the torsion would lead to the following equation, instead 
of (3.12): 

(1) „ .(1) (1) 



{<t^r - 0bg) + T. (3.23) 

Here T is a suitable density caused by the non-triviality of the torsion components, say, 
= —T21 (this is just appropriate for the screw dislocation along Oz^ [51]). The notation 

(1) 

(jy-j- implies the corresponding solution at a given T. Being considered as an additional 
source of incompatibility. T can be assigned to possess near Oz a series expansion with 

(1) 

arbitrary, though adjustable, coefficients. As a result, ^-j- can also acquire, in comparison 

(1) 

with (3.13), certain modifications near p = 0. 

Therefore, (A 0j is replaced, on an extended treatment, by (A (pq- ) . This can be 
effectively accounted for in the present framework as well. Let us introduce a piece-wisely 

constant density to "regularize" ( A j as follows [50] : 



where /i[o,p,](p) is unity for p e [0, p*] (i.e., within a disc) or zero otherwise. It is crucial that 
(3.24) should not be treated as a replacement to be iteratively improved. Let us use (3.24) 
in (3.16.1). Then, an analogue of (3.19) is valid, where the corresponding "potential" ^ is 
present: 

9 = ^(|-i)Wp)+C- (3.25) 

b b 

The constants C and p* should be specified later. Both —k^Kq{kp) and — -h\o 0,1 (p) are 

27r npi 

5-like for k and 1/p* large enough, respectively. Considered as the surface densities, they 
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are properly normalized. Therefore, the replacement (3.24) looks better, provided that its 
S-\ike character is sharper. 

Using (3.19), where given by (3.25) is substituted, from (3.15) we obtain the following 
equations [50]: 

(A - k') (a + (/-/bg) = kn, (3.26.1) 

P = -T^ A/ - (M/33 + 9) + T^{f- /bg) . (3.26.2) 

I — a 1 — a I — a 

Here, /bg respects the first equation (3.20), and k is given by (3.21). The driving source in 
(3.26.1) is defined as follows: 



7^ ^ (A-K^){^,,+g-^^,l-Q') 

(3.27) 

+ (l-a)(l - 2f,'Cr){^ - $bg) - ^(A 



It is assumed that \1'33 and $ (or, analogously, \E'33 and $bg) in (3.27) are given by (3.18) 

(1) (1) (1) 2 

where (or, correspondingly, 0bg) is substituted appropriately. In addition, (A 0) is 

kept in (3.27) only formally, as a symbol. The first equation in (3.26.1) (would be called, 
after [29] and [30], the non-homogeneous bi-Helmholtz equation^ is just to be solved in what 
follows to determine the stress function /. Equation (3.26.2) is needed only to express the 
stress component azz = P- 

The governing equations (3.26) are essentially similar to equations that should be ex- 
pected under a consideration extended by means of the torsion (see also explanations in 
Section 5.3). The actual derivation of (3.26) incorporates (3.16) and (3.17), which are sub- 
ject to (3.24). Therefore, the teleparallel framework is valid only approximately. However, 
the free parameter is coming to play by means of (3.24), and this supports the simi- 
larity noted. Indeed, the arbitrariness owing to is reminiscent of a freedom caused by 
the non-triviality of appropriate torsion. Technical complexities, which might appear in the 
description based on (3.23), should be avoided now, since they deserve separate studying. 
Therefore, the strategy based on (3.24) looks appropriate. Equations (3.26) lead eventually 
to a self-consistent stress distribution without an artificial cut-off at the core. 



3.3 The driving source of the gauge equation 

Before solving (3.26.1), simplifications for TZ given by (3.27) are in order. It is convenient to 
investigate TZ multiplied by p^. We use (3.13) and (3.18) and keep the same notation TZ after 
the multiplication. After re-scaling the radial coordinate up t— > s, we obtain approximately 
(see [50]): 

3 

7^ ^ s^w{s) + ^7^a, (3.28) 

a=l 

^ See [29] and [30] for appropriate refs. and more examples of non-homogeneous bi-Helmholtz equations. 
A 5-driven differential equation analogous to (3.26.1) has already appeared in [15] as well. 
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where 



X-27^2 ^ -csK,{s)As(l> '^sK,{s)h[o,s.]{s), (3.29) 

d 



The following notation is accepted in (3.28) and (3.29): 

c = (l-a)(l-2/x2C7) - 4c, c = iJ,^C3, 



2 r I 



(3.30) 



where a and 1 — a are given by (3.10), s* = /tp*, and /i[o,s,](s) is unity at s G [0, s*] or zero 

(1) 

otherwise (we use C = Q' in (3.25)). Further, in (3.29) implies with removed multiple 
b /6k\2 

— . The factor X'^ = I — ) just points out the fact that the driving source is quadratic in 

d d 

b. Moreover, is s s ^ ), and Kq and Ki are the modified Bessel functions [59]. 

ds\ ds/ 

(1) 

The following comments on (3.28) and (3.29) are in order. Classically, A 0bg is the 
Dirac's delta-function, and its square would occur instead of (3.22). The equation for the 
corresponding stress function /bg is conventionally considered for strictly positive p. Then, 

(1) 

the corresponding driving source is not influenced by A ^^g. Therefore, it seems inappro- 
priate to account for, at the same footing with (3.24), the contributions expressible in TZ by 

(1) (1) (1) (1) 

the differences A —A ^p(^ 0bg)- To put it differently, the regularization 

allowed classically "persists" into the non- conventional approach, and (3.28) and (3.29) take 

their actual form. 

The asymptotic properties of TZ (3.28) are as follows. First, TZ is well localized since it 
contains the modified Bessel functions. At small s, we expand 

7l{s) ~ pi + P2 log S + P3 

+ (p4 + P5 logs + pe) + o{s^) , 

where the numerical coefficients pi,p2, ■ ■ ■ ,P6 are 
Pi = -X'^c, P2 = 0, 



(3.31) 



P6 = 



«,(0) + X^(l-21ogI)((l-21og2-l)± + (c<-l)|), 



(3.32) 



with a = 1 or a = 3 in p3,p5, and pq. Indeed, the calculation of 7^3 requires, as soon as (3.24) 
is imposed, the differentiation of the step-function. This would result in singularities in the 
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driving source. Our regularized scheme can thus be destroyed. We use P3,P5, and pe with 
a = 1, when TZs is simply omitted to express the neglect of the corresponding differentiation. 
Otherwise (for comparison), a = 3 when TZ^ is calculated according to (3.29). However, a is 
left unspecified to give an indication of both possibilities. Moreover, in spite of p2 = 0, it is 
instructive to keep the corresponding term in (3.31). It is to be equated to zero at the end. 



4 Solution of the gauge equation 



We are going to solve (3.26.1) in two steps [50]. First, we consider the non-homogeneous 
Bessel equation 

d 



ds"^ ds 



G(s) = 7^(s) . 



(4.1) 



The method of variation of parameters [62] provides a general solution to (4.1). We choose 
the modified Bessel functions /q and Kq [59] as two hnearly independent solutions of the 
associated homogeneous Bessel equation. Solution to (4.1) appears as follows: 



G{s) = lim 



where 



s 

Io{s) (A(e) + J Komit) J 

€ 

S 

+ Ms) (5(6) - J lomit) j) 

€ 

oo 

A{e)^ -jKomit)j, 

e 

const + |(| + |logt+(^+p3)^) 



(4.2.1) 



(4.2.2) 



dt. 



The behaviour of t~^/o(t)7^(t) and t^^KQ(t)TZ(t) for small or large t is important for 
justification of (4.2). When t is small, the regularization in G{s) is ensured by means of 
the specially constructed A(€) and S(e) from (4.2.2) (see Appendix A). When t is large, 
Ko{t)TZ{t) decays as e~^*, whereas Io(t)7l(t) behaves, mainly, as a constant. Therefore, the 
infinity as an upper bound in A{e) given by (4.2.2) is appropriate. The choice of the integra- 
tion bounds implies that G{s) decays exponentially for large s. The important arbitrariness 
in B{e) from (4.2.2) should be stressed, which is due to the additive constant denoted as 
const. 

Using Appendix A, it is straightforward to expand G{s) for small s: 



G{s) ~ gi s 



-2 



3 3 

+ Yl log' s + s'^Y^ qg^i log' S + 

i=0 1=0 



(4.3) 
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The numerical coefficients in (4.3) are expressed in terms of the parameters (3.32) as follows: 

Pi P2 Pi . P3 

3»i , 7 _ pi , . . , 

?4 = - const + ^ > ?5 = - const x log - - - — , (4.4) 

P2 Pi , P3-P2 , P4 

^^-24' ^^-32 + ~8~ + T' 

where Ik is given by (A8) in Appendix A, and const in q4 and is introduced by (4.2.2). 
The coefficients gg and qg in (4.3) are practically too complicated and are not of importance 
below. 

As a second step, we express the difference / — /bg = k J-', where J-' respects the non- 
homogeneous Bessel equation 



ds'^ ds 



-.2 



The modified stress potential of the second order appears as / = /bg + k T. For convenience, 

(4.5) is also written in terms of the new variable s, which is however defined differently: 

^2 

s = J\fp, where N"^ = . We choose Yq and Jq as the fundamental solutions of the 

1 — 2a 

homogeneous Bessel equation associated with (4.5) and obtain 
^(p) = CYoiAfp) + DUNp) + /^(p), 

I^{p) = Jo{Np) J Yo{Nt)G{Kt)tdt - Yq{Mp) J JoiXt)G{Kt)tdt, 
p p 

where G{-) is given by (4.2), and Yo{s) = {7r/2)Yo{s). Here the integrals are convergent at 
their upper bounds. 

In the preceding consideration (see [50]), the choice of C 7^ and D = has been made in 
J^{p) from (4.6). However, in this case the stress field looks non-conventional for sufficiently 
large distances. In the present paper, we put C and D equal to zero, and thus J-'{p) = Ir{p)- 
The integral Ij^{p) decays exponentially for large p, while for small Hp we obtain (Appendix 
B) 



T{p) ~ ro + ri logp + p^ X^'^e-i logV -h ryp^log^p . . . , (4.7) 

where 



1=0 



ro = -Xy + log^A/") Xj , n ^Xj, r2 = ^ 

P2 (s Pi /1 , ^P2 . P3 

1^3 = w::, U 



24 ' 



/ 2 \r2\ P2 
''^^^ --^^384' 

and pi,P2, and pa are given by (3.32). The coefficients and re are given by (B6)-(B11), 
and Xy and Ij are defined by (BIO) in Appendix B. 
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5 The stress tensor 



5.1 The components Gpp and cr^^ 

Therefore the modified stress potential of second order / takes the form 

f^h^ + ky, (5.1.1) 

fbg ^ -n log^ p + di + (^2 logp , n = k — , (5.1.2) 

where /jf is given by (4.6), and the stress potential /bg is in full agreement with [7]. The 
potential /bg respects the following equation, which results from (3.20) for p ^ 0: 

AA/bg = ^ . (5.2) 

Taking into account (3.8) and (3.11), we represent the stress tensor of the second order 
in the cylindrical coordinates as follows: 



(2) Id (2) £ (2) 



(2) 

where / from (5.1) and p from (3.26.2) are substituted. The other components of a are 
zero. In the classical problem, we use (5.3) with /bg from (5.1.2) and obtain 

' ' (5.4) 
.(2) o~logp , 2n + o?2 ^, 
(<7bg)</.0 = -2n — 5 2di. 

(2) 

The free parameters di and ^2 in (5.4) are given by the vanishing of ( a hg) pp on the boundaries 
of the hollow cylinder, e.g., p — Pc and p — Pe > Pc'- 

{'^hg)pp = hg) pp = 0. (5.5) 

p=Pe p=Pc 

Substitution of (5.4) into (5.5) allows us to define di and ^2- For our purposes, it is 
more appropriate to assume that Pe is essentially greater than p^ {pe / p is moderate). We 
approximately obtain 

di ^ log — , (i2 ~ 2nlogpc , (5.6) 

Pc pc 

and _ 

.(2) . 2nf p fP\\Po\ 
{(^hg)pp ~ — log ( — ) log— , 

P^ V Pc ^Pe ^ PcJ ,^ 

~ (5.7) 
.(2) . _ 2n/ p fPY. Pc^ 
f^bg 00 ~ — 1 - log — ] log — 

P^ V pc ^Pe ) pc 

Note that n is proportional to 6^, and /bg from (5.1.2) is indeed of second order in h. From 
(3.29), (4.2), and (4.6) it is seen that Ijr is also quadratic in h. 

Let us turn to the stress components (5.3) with / substituted from (5.1). The solution in 
question depends on several free parameters: p*, const, di, and ^2 (see (3.25), (4.2.2), and 
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(2) 

(5.1.2), respectively). Moreover, a from (3.26.2) depends on C (= Q'). In order to avoid 
the arbitrariness, the asymptotical stresses should be subject to appropriate conditions at 
rep ^ 1 and Kp <^ 1. For instance, Ij: is exponentially small for «;p ^ 1, and therefore 

/ ~ /bg(p)- As a result, ^(7 pp— (^(^hgjpp arid ^(t {^o^bg)(t><t>^ where the background stresses 
are given by (5.4). However, di and o?2 in the gauge case are defined below. 

(2) (2) 

In the limit Kp <^ 1, we use (4.7) and (5.1) and obtain a pp and a as follows: 

app-{%\,)ppc^ - 2r2i^-^-2r3logV-(3r3 + 2r4)logV 

- 2(r4 + rs) logp - (rg + 2r6) - ^ttp^ log^ p + . . . , 
(2) .(2) . „ logp 2r2 - ri (5.8) 



^2 „2 



p- p- 

- 2r3log^p - (9r3 + 2r4) logV - 2(3r3 + 3r4 + r5)logp 

- (3r5 + 2r4 + ^tq) - Urrp'^ log^ P + ■■■ , 

where the background stresses on the left-hand side are given by (5.4). Now, ri, r2, 
and Tj include, for notational compactness, the multiple k from (5.1.1). The contributions 
oc p^^logp disappear on both sides of (5.8) since r2 is equal to n. The contributions due 
to r3 and disappear in (5.8) since p2 = in (3.31). However, the terms, which are either 
divergent or constant as p — > 0, are still contained in (5.8). 

The number of the free parameters is smaller in comparison with [50] since now C — 
in (4.6). Therefore, the inappropriate terms in (5.8) should be handled differently. Recall 
that the core region corresponds to p < k^^. Let us introduce a fictitious boundary at p = p^, 
where p^ -C We assume that p varies within the segment [pe,Pe], where p^ is allowed 
to be arbitrarily small but strictly non-zero. Then we arrive at the following equations: 



ri + d2 



+ 2r4 log' p, + 2(r4 + rg) logp^ + rg + 2(r6 + di) = , 



- 2r4 log' p, - 2(3 r4 + r^) log p, - 3 rs - 2(r4 + rg + di) = , ^^"^^ 

Pe 

2nlogpe - (i2 - 2p2 (ii = 0. 
The first two equations in (5.9) imply that the terms, which are either divergent or constant 

(2) (2) 

in (5.8), vanish at p = p^. Therefore, a pp and a at p — p^ are mainly given by the terms 
that go to zero as p^ — > 0. The third equation in (5.9) reads that the radial stress is zero at 
the outer boundary p = pc . The arbitrariness of di, ^2, and const can be used to meet (5.9). 
Since di, ^2, and const are thus determined (as functions of p^, in fact), we can use them 
in (5.1) and (5.3) to produce the stress distribution valid on the ring, say, I^[p£,pe] confined 
between p — Pe and p — Pe ■ 

Let us consider an infinite sequence of concentric circles with monotonically decreasing 
radii p = p^ (pe's remain however non-zero). Each of the corresponding rings ^^[p^,^^] is 
completely covered by a disc with the same external boundary p = Pc but with a puncture 
at p = 0. The stress distributions defined on the rings also form a sequence. It can naturally 
be considered as converging to an appropriate solution defined on the punctured disc. Such 
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go to zero for decreasing 



(2) 

a limit is expected to exist, since the boundary values a.^ 

Pe- Therefore, our attention is focussed at the rings possessing sufficiently small internal 
radii. Note that even the definition (4.2) is concerned with a hollow domain possessing a 
sufficiently small internal radius. It is just the regularization in (4.2), which allows us to 
send this radius to zero immediately. 

To extract the dependence on const in (5.9), let us turn to ^4 and ^5 from (4.4). Using 
Appendix B, one can establish that ri in (4.8) and r5, re in (Bll) depend on const as follows: 

2 2 
ri = -C + Fi, r5 = -^C + ?5, r6 = ^C + f6- (5.10) 



The following notation is accepted in (5.10): 

~ k 
C = const X 



k2 + j\^2 ' 



log(|Ac), (5.11) 



where k'^ — /i/l {/i is the shear modulus, and i arises in (2.10)) and {k/MY — 1 — 2a. 
Moreover, a and k are given by (3.10) and (3.21), respectively. We substitute (5.10) into 
(5.9) and rewrite the result as a single 3x3 matrix equation: 



(5.12) 





where 



■1 + 



4 



(5.13) 



2 

h = 2ripl\ogp^ + {ri + n)pl - 
I2 = -2ripl\og^ pe - 2(r4 + r5)p^logp^ - (F5 + 2r6)p^ - r 1 , 
/s = 2nlogpe ■ 

Equation (5.12) is solved for di, ^2, and C by the Cramer's rule: 

c = I {pI {h - k) + plik - h)) , 

di = ^ {ak - bk + h{b - a)) , (5.14) 
2 

c^2 = ^ (Pe (^^1 - ak) + plak) , 

where V = 2a pi + 2(6 — a)p^ . The parameters C (i.e., const, see (5.11)), di, and d2 thus 
obtained should be used in / from (5.1). 

Let us assume that p^ is sufficiently small (the so-called logarithmic approximation) : 

11,1 /X 

1<K< — <— log— . (5.15 

Pe Pe Pe 
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Then, the ratios (5.14) can be represented in the perturbative form 

C ~ -^log - + 4(^5 + (1 + 5)^4) + ... , 

~ + 2^1og -)+... , (5.16) 

Pi V Pi pe' 

^4 1 _ 1 

(^2 ^ - 4— log 2n log h . . . , 



where dots imply terms that tend to zero as — > 0. In addition, the new length parameter 
Pi is intended to express the contributions which are constant as — > 0: 

~ 1 ~ 4 \ 
2nlog— = ri ^r5 + (l + g)r4 , (5.17) 

P\ / 

where q is given by (5.11). Using (5.4) and (5.16), we obtain sufficiently far from the core: 



(logii-(iLfiogM_4lilogpii-4), 

V Pi Vpe/ Pi/ Vp^ Pe^ 

- log^ - (^)^og^) + 4riiogp.(l + . 

p^ V Pi ^Pe ^ Pi/ «^ ^P^ Pe ^ 



(2) 2n 

^ ^ ''^ ' ''^ (5.18) 

(2) ' . . . . 



As far as the asymptotic expansions (5.18) are concerned, it is crucial that r4 can be 
made equal to zero. Indeed, using in the form (4.8) (where pi and p3 from (3.32) are used, 
and p2 is zero), we obtain p^ from the equation r4 = as follows: 



2_5(^,") 8(l-i.)(l-r^)_^ 

(5.19) 



p* = — zr^ — ' ^(^'«) = 



K'- q;(1 — i/)?7 — 2i/(l — 77) ' 



?7 = 4(1 + v) 



2p2C7 ' 



where a = 1 or a = 3, according to (3.32). Equations (5.19) lead to the positivity of the 
coefficient-function B{r],a). The vahdity of (3.24) holds better whenever k,~^ and p* are 
smaller. If the scale is small enough, a bound Bq can be found to restrict B{r), a) : 
< B{r],a)^ Bq. If so, p* from (5.19) is estimated as p* = 0{k,~^). In terms of 77, the 
corresponding restriction takes the form 

where < 1/ < 1/2 for isotropic materials. 

It is straightforward to rewrite (5.20) in terms of the elastic parameters C3 and C7: 

a) Cr + 2{l + u)C3 <^<Cr + 2{l + u){l + U)C3, C3 > , 

" ^ (5.21) 

b) C7 + 2 (1 + z/)(l + [/) C3 < — < C7 + 2(1 + u)C3, Ca < . 

Further, (5.21) can be re-expressed [7] in terms of the parameters m' and n' from (3.4), which 
can in turn be related [4] to the third order elastic moduh of the cubic crystals. In other 
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words, the choice of material is restricted by (5.21). The numerical data provided by [4] 
witness that these conditions look reahstic for certain materials. However, to discuss (5.21) 
from the viewpoint of real crystallography is beyond the scope of the present investigation. 

Straightforward usage of (4.6), (5.1), and (5.3) leads us to the following general expression 
for the stress components: 



(2) _ .(2) . , _iw 

pp — I cr farjpp — ftp ijp , 

(2) (2) 

= (o-far)<^0 + k{M^Ir + p~^I'r - G) , 



(5.22) 



where k and G are given by (3.21) and (4.2), correspondingly, and I'-p (the prime implies the 
differentiation d/dp) is expressed as follows: 

oo oo 

N-^I'^ = Yi{Np) J Jo{Nt)G{Kt)tdt - Ji{Np) J Yo{J\ft)G{nt)tdt. 
p p 

It is assumed that Yi(-) = (7r/2)yi(-) (similar to Yq in (4.6)). Now, the subscript 'far' labels 
the stress components in (5.22) in order to emphasize the non-conventional choice of di and 
(i2, which corresponds to (5.16), provided that the internal radius of the ring is small 

enough. We also obtain the sum of the second order stress components: 

(2) (2) (2) 

+ ^H= (^far)pp + (C^far)^<^ + ^(AA^/^ - G) . (5.23) 

Since r4 is zero for pi from (5.19), the limit ^ is suitable in (5.16) to eliminate 
Pe. Then, di and d2 acquire the classically known form (5.6), though pi is present instead 
of pc . The limiting value of the parameter const (see (4.2)) is also expressed with the 
help of (5.11) and (5.16), provided that p^ decreases to zero. The large-distance stresses 
(5.18) take at r4 = the classical form (5.7). Note that the axis Oz is not captured by 
the present approach, as far as the limiting transition is allowed. This is in agreement with 

(2) 

the regularization of the density profile at p = in Section 3. The components {(Tiax)pp 

(2) 

and (o'faiOtiij reduce, after the shrinking of the internal boundary of the ring, to (5.18) (with 
r4 equated to zero) and should be distinguished from the analogous ones given by (5.7). 
Eventually, there is no cut-off around the dislocation's core for the second order stresses 
expressed by (5.22), provided that the numerical parameters are substituted appropriately. 



(2) 
0-pp 



5.2 The component Gzz 

(2) 

Let us consider the stress component (J zz— P- Classically, the second equation in (3.20) 
implies the following representation [7]: 

(2) (2) (2) 

<y zz = y (<ypp + 'Jh) - 2/^(1 + ^) (^'f + Q') > (5-24) 

where /bg from (5.1.2) is used for the stresses on the right-hand side (the respective subscript 

(1) 

is omitted for brevity), and 0bg= {~^/'^^) ^ogP is substituted to express ^^^f. As far as the 
equation (5.24) is concerned, it is assumed that the argument p varies within the segment 
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(2) 



[pcPe], and Q' is determined from the requirement that cr^^ averaged over I^[p^,pg] is zero 
(the so-called 'mean stress theorem'; see [4]). Then, the representation (5.24) reads for 
Pe > Pc : 



(2) 
(7 , 



log — 



Pi 



(2) 



Let us turn to the stress component o from (3.26.2) of the modified defect: 



(2) 



(2) (2) 

y (V,, + a^^) - 2/x (1 + !/)(*,, + 5) + (1 + (/ _ 



(5.25) 



(5.26) 



(2) (2) 

Now p varies within [p£,pe]- From the previous considerations it is clear that + <^,^ 

(2) 

and ^22, taken at p = p^, decrease as 0. Nevertheless, a behaves, sufficiently close 

to Oz, somewhat artificially. The reason is that the logarithmic terms dominate in / — /bg at 
p small enough. However, this is not an obstacle to view the rings 'D\p^^p^\^ as well as (in the 
limiting sense) the punctured disc, as the domains of definition of the stress distributions. 

(2) 

Indeed, o from (5.26), averaged over I^[ps,pe], is finite, provided that p^ goes to zero. 

(2) 

Let us obtain a ^z from (5.26) at sufficiently large distances, where / — /bg (= /c/jf) is 
exponentially small. Using (5.18) at r4 = 0, we obtain 



(2) 



cr 22 ~ 2i/n 



1 2 pe 

— - — log — 

P Pe Pi. 



- 2p(l + v) 



We determine Q' approximately: 



Q' 



( b 



log(Pe /P2) 



Pi 



(5.27) 



(5.28) 



where another length p2 is introduced, namely, 

00 



6 \2 C 



(5.29) 



(; 



pdp + 2k^ — — — / Ij^ pdp , 
l-v J 

1 
0, and the upper integration bound pe is replaced approximately by 



provided that pe 
infinity. 

Therefore, (5.27) depends on two scales, which are, in principle, different in our consid- 
eration: pi from (5.17) and p2 from (5.29). To re-write (5.27) more conventionally, it is 
appropriate to introduce another length, say pm , as follows: 



vn 



/ 6 \2 



log 



_ 1 f b 1 

= vn\og — -\- — iJiil-\- v)c\og — 
Pi V27r/ P2 



(5.30) 



pi \27r. 

Substituting Q' from (5.28) into (5.27) and using (5.30), we obtain 



(2) 
(7 , 



yn + i^—j /x(l + i/)c 



P" 



2 , Pe 

log — 

Pir 



Pe^ 



(5.31) 
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Equation (5.31) looks similarly to (5.25), except that now pm is present instead of pc • 

The classical answers for di and 62 are given at pe ^ Pc by (5.6). Since di and 62 from 
(5.16) include the constant contributions, the agreement between (5.6) and (5.16) is due to 
the newly defined length pi from (5.17). The mean stress theorem implies that p2 coincides, 
classically, with the lower integration boundary pc , provided that pc remains unspecified. 
Instead, Q' in (5.28) is concerned with another length p2 from (5.29). This is because zz 
is unconventional inside the core. In the present approach, two different lengths arise. The 
definition (5.30) is only to obtain the classically-looking representation (5.31). 

We use 'g (3.25), where C 
expressed by (5.26): 



(2) 



Q' with Q' given by (5.28), and we rewrite the stress a 
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.(2) 



cTpp + - 2p(l + z/) 
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7: ^ 



Pi 



m^z p dp 



2p 



+ 2k' 



(5.32) 
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Moreover, using (5.23) we can further rewrite (5.32) for sufficiently large pe as follows: 

4) - 2^(l + z/)^,, 
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- 2p(l + z/)( 
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cr far) 
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log 
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(5.33) 



+ k [M^y - vG), 
where fc, G, Jjr, and pi are given by (3.21) 



for ^! zz from (3.18) in the polar coordinates takes the form: ^! ^ 



(4.2), (4.6), and (5.19), respectively. Expression 

(1) „ (1) 



-c{dp 0)^, where is 



given by (3.13). The behaviour of ^ zz for sufficiently large distances is similar to that of \Ef^f 
although is different within the core. 

Let us recall that a stationary Schrodinger equation obtained in the effective mass approx- 
imation [63] has been used for the conduction electrons in the dislocated [13] or disclinated 
crystals [44], [45]. The influence of the defects has been accounted for via the deformation 
potential given by the trace of the strain tensor (i.e., by the dilatation). The gauge and/or 
geometric approaches to various effects due to the dislocations (due to the torsion) attract 
considerable attention: [64]- [70], [72], and [730. 

The dilatation is non-trivial for the screw dislocation just because of the strains of the 
second order. Sufficiently far from the core, the dilatation depends on 1/p^ just by means of 
m' and n' from (3.4) (see [7]). The results obtained above can be applied to the Schrodinger 
equation, as is discussed in [13]. Indeed, the solution [13] is valid only outside the core, and, 
say, a traction-free condition is still required. In turn, a possible influence of the dilatation, 
obtainable by the present approach, on a local shape of the wave function can be investigated 
in the vicinity of the core region. The corresponding results should be comparable with 
appropriate lattice-based investigations. Note that the disclination core's radii have been 
used in [44] and [71] for description of the corresponding electron localization. 



*For more refs. on other types of defects, see [57]. 
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5.3 Remarks on the logarithmic approximation 

The asymptotic relations (5.16) give an indication of the fact that the ratios (5.14) should 
be re-arranged perturbatively for log(l/p£) and pe sufficiently large. It is instructive to re- 
derive the leading terms in (5.16) by means of a more straightforward use of the logarithmic 
approximation at r4 = 0. Indeed, let us reconsider (5.12), provided that the parameters 
(5.13) are taken in the principal order as follows: 

a fti -1 , b ^ -1 — logpe , 

2 (5.34) 

h ^ -ri, k ^ -ri - 2 r^pl log . 

Here, r4 = in ri, and the special notation is chosen for at = 0. Besides, is kept 
like in (5.13). We replace in (5.12) the entry equal to 2p^ by zero. Then, in agreement with 
(5.16), the solution reads: 

4 ° ~ 

C = , di = log — , d2 = 2nlogpi, (5.35) 

Pi Pi 

where pi agrees with (5.17). 

Moreover, one can estimate pi, provided that logK is large enough. With regard to the 
definition by means of the equation (5.10), we extract the dependence of ri and re on the 
large logarithm as follows {r^ is zero): 

ri = log K + ri , tq = r5 log k + tq. (5.36) 

The terms ri and tq in (5.36) are simply intended to denote the corresponding remnants, 
which are finite, provided that log k is growing. We can use (5.36) in (5.13) and again directly 
solve (5.12). Then, di and ^2 appear just in the form (5.35) with pi — Bi/k, where 

2nlogSi = -ri + ^rs. (5.37) 

It is seen that (5.37) agrees with (5.17). The radius p2 can also be estimated. Prom (5.29) 
we obtain p2 — B2/K, where 

1 

2clogS2 = - \ -2c /(I - sKi{s)y—, (5.38) 
47 s 


and c is given by (3.30). Therefore, the logarithmic approximation allows us to determine 

the leading behaviour of C, di, and d2, as well as to estimate pi and p2- Note in passing 

that a certain infiuence of the regularization in (4.2) on Bi from (5.37) and B2 from (5.38) 

should be accounted for, though elsewhere (since specific numerical estimates based on the 

crystallographic data are required). 

The result (5.35) should be commented as follows. Let us introduce a sufficiently small 

positive G by the requirement: absolute values of the stresses, which are smaller than G (say, 
(2) G, 

I CTij I < 77), should be treated, for engineering, as indistinguishable from zero. Then, the 



23 



(2) 

shrinking of the boundary at p = can be ceased just for those boundary values a ij 

P=Pe 

which respect this estimate. At the same time, (5.15) holds, and thus (5.35), (5.37), and 
(5.38) are valid. The large distance behavior looks conventional though includes now pi 
and p2, which are dictated by a choice of material. Provided is thus fixed, the stress 

(2) 

a zz remains bounded outside an infinitesimally thin tube around Oz. It is essential that 
under the prescribed accuracy e, a distant "observer" outside the core is dealt only with the 
physical stresses improved by means of pi and p2- These stresses are insensitive to particular 

Let us turn again to our geometric interpretation. Let us imagine for a moment that the 
differential-geometric torsion is taken into account in order to make (3.16) conventionally 

(1) (1) 
looking. Then, A ^ comes to play instead of A . Since a specific shape of appropriately 

(1) 

localized density T is unclear, (3.24) can be applied, as a simplification, just to A 0^. As 
a result, the driving source TZ should be changed. The local properties of this new TZ can 
be estimated. The modifications are expectable for p4, ps, and from (3.31). However, 
these changes do not influence the coefficients we consider in (4.7). A change of shapes of 
the profiles of the stress components is expectable in the middle of the core, but this is 
beyond our scope. It is important that pi in (5.17) and p* in (5.19) should not be influenced 
(1) 

when A 0-7- is subject to (3.24). Only p2 from (5.29) should be replaced, since it is deflned 
"globally" . Thus, the present picture is almost the same as that with the torsion admitted 
but handled approximately. To avoid a specification of T, it is appropriate to keep our 
teleparallel interpretation. 



6 Discussion 

The stress potential of the non-singular screw dislocation is found in the quadratic approx- 
imation. It is given by the sum of the conventional part and the gauge contribution. The 

(2) (2) 

general expression (5.22) for the stress components of second order a pp and o",^,^ is also 
obtained by means of appropriate differentiations of the stress potential. In addition, the 

(2) 

general expression for the stress o is elaborated as well: see, for instance, (5.32) or (5.33). 
The main attention is paid to the asymptotic properties of the stresses in question. To ensure 
their self-consistency, the arbitrariness of certain constants in the general solution for the 
stress potential is used. 

(2) (2) (2) 

Sufficiently far from the core, the second order stresses obtained, a pp, a and a ^z 
are in agreement with [7], although now they include the new lengths pi and p2. The non- 
conventional part of the total stress potential is responsible for the localized contribution 
into the full stress distribution. As a result, the modification of the short-distance behaviour 

(2) (2) (2) (2) 

is as foHows: a pp and a tend to zero, while a ^z grows logarithmically (although a is 
integrable, as required in [7]) in the close vicinity of Oz. There is no artificial cut-off at an 
internal boundary. The "exterior" (with respect to the core region) solution for the second- 
order stresses plays a central role in the picture presented. This is because of its improvement 
due to the self-consistent arising of pi and p2. The general solution is characterized by one 
more length p^, which is a half-width of the effective defect's density profile. 

The lengths pi, p2, and p* are expressed through the elastic moduli of second and third 
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order. The appearance of several lengths correlates with the absence of a sharp boundary 
around the core in the first order pictures in [14], [15], [17], and [19]. In other words, a 
concept of 'transition shell' looks helpful for characterization of the solution obtained. This 
shell would separate two domains: a compact region under the shell, where the strains are 
rather finite, and an outer bulk. The solution obtained should be less reliable under the shell 
(where the effects of discreteness are comparable with those of strong elastic non-hnearity; 
see [3] and [4]). Within the transition shell, the quadratic corrections both of the conventional 
and of the gauge origin (i.e., due to Ij: in the stress potential) should be important. Outside 
the shell (i.e., in the bulk), the conventional terms become more valuable, since the gauge 
contributions are strongly decaying. The radii obtained pi, p2, and should characterize 
the location and the extent of the transition shell. More generally, they seem to display the 
dislocation cores as radially-layered regions. 

The radii p*, pi, and p2 include the basic length 1/k { = ^ which can be adjusted 
to the interatomic scale even in the first order: see, for instance, [14] and [17]. Analogous 
identifications can also be found in [20], [24], and [26]. Here ^ is the gauge- material parameter 
(see (2.10)) and p is the shear modulus. However, the corresponding coefficients B{r],a), 
Bi, and B2 given by (5.19), (5.37), and (5.38), respectively, are strictly influenced by the 
crystallographic moduli of second and third order. MutTial comparisons of the coefficients and 
more specific statements on the properties of the transition shell should be done elsewhere. 
The numerically obtainable values for pi, p2, and p* can be compared, in principle, with 
appropriate scales provided by the semi-discrete models of the cores (see refs. in [3], [4], and 
[10]). On the other hand, attempts at matching the radii in accordance with the experimental 
observations may result in independent estimate for Such estimate can, in turn, be 

compared with 1/k obtainable by means of the ffist-order considerations. 

As far as the numerical estimates are concerned, the isotropy requirement should be taken 
into account, since the number of the elastic constants of third order for the cubic crystals 
is six (see [4] and [10]). Moreover, strong couplings could modify the elastic parameters 
within the core in comparison with those outside it (see [40]). Note that [36] also provides 
the radius of the screw dislocation as a function of the elastic constants (of second order 
though) and of the lattice spacings of the crystal. 

Strains and displacements within the transition shell can, in principle, be elaborated with 
the help of the present approach. For instance, the displacements of the first order in the core 
region are discussed in [14] and [17]. The next order corrections to those displacements can 
be obtained within the shell. Here, the technical issue [74] could be helpful. Thus corrected 
displacements should be comparable with the semi-discrete calculations subject to the so- 
called flexible boundary conditions (see [75] and [76]). Direct observations of the dislocation 
cores by means of the high resolution electron microscopy can also be used for comparison 
(for more refs. on the observation and modelling of the dislocation cores at atomic level 
in semiconductors, for instance, see [78]). The solution obtained can also be tested for the 
problems where just the second-order elasticity is relevant, e.g., for the elastic waves or the 
electron (see Sec. 5.2) scattering, etc. (see [4] and [9]). 

The continuation obtained for the stresses is just due to the gauge equation considered 
as the incompatibility law. In turn, our constitutive equations imply that the elastic energy 
is written up to the third order. But the quadratic constitutive law is only an approximation 
for essentially non-linear situation. Clearly, the present picture does not pretend to simulate 
the stress distribution in the middle of the core, where both the atomistic structure and 
the finite elasticity are crucial (see [77] where a non-linear elasticity approach with fiexible 
boundary is discussed as a way for improving the semi-discrete techniques). Nevertheless, the 
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stress distribution provided by the present model of the screw dislocation is self-contained 
outside a certain shell, and it seemingly remains satisfactory within this shell as well. 

Our investigation demonstrates that the Hilbert-Einsteinian gauge approach is flexible 
enough, since it allows, in two orders, the self-consistent description for the non-singular 
screw dislocation. Certain features of an extended picture could be discovered in the present 
model. It can be used as a base for further, more involved, treatments. Application to the 
edge dislocation would be also desirable. 
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Appendix A 

Appendices A and B contain intermediate results, which are helpful for expanding the stress 
potential / in (5.1) (see also [50]). Starting with the series expansion (3.31) and using the 
Bessel functions in the series form [59], we go over to the corresponding series T from (4.7). 

First, let us estimate G'(s) from (4.2). To begin with, we expand the products 
t'^ Ko{t)TZ{t) and Io{t)7l{t) for a small argument t <^ 1: 

rl/o(^)7^(^) ~ pir^ + pst'Mogt + + ps) 

(Al) 

+ (p4log't + ( J + Ps) logt + ky + ...- 

t-^Ko{t)n{t) ~ -Pi \og(^t) t-^ - (pslog't + ki \ogt + k2) t-^ 

^2 / ' (A2) 

- (p4 log^ t + ks log^ t + kilogt + k5)t + 

where Pi,P2, ■ ■ ■ ,P5 are given by (3.32), and the coefficients ki, k2, and k^ are expressed as 
follows: 

^1 = ^ + logQ)p2 + P3, 

-7 - '°^(i) (t - • («) 

ks = — + log(^-jp4 + P5 

(7 is the Euler constant). The parameters k^, k^, and k are not used in the present paper 
and can be found in [50]. 

^The author is grateful to M. Lazar for calling attention to [79] where the modified screw dislocation has 
also been obtained for the gauge Lagrangian of the Maxwell type. 
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Using (Al) and (A2), we obtain the following estimates for s <S 1: 



s 



(/Ci logs + JC2) 



+ ^ /Ce-i log* s + ^ /Cio-i log* s + ... , 



i=0 



where 
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+ P3, 



= -constxlog^ + (1 + 31og^)^, 

- (1 - log 
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+ log( 
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+ log's + ^Jio-i log*s + 
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i=0 



where 
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-Pi 



2 ' 
-Pi 



8 



+ k2, 



^ P4 

12 ^ 2 ' 



'P2 

3 ' 

Xe=X.+ (l.21„gI)f, 



(AT) 



and ki, k2, and k^ are given by (A3). The numerical constant Ik in ^6 from (A7) implies 
the regularized value of the integral: 



Ik = / Ko{t)n{t) 



dt 



Ko{t)n{t) + Pi iog(|t)r 



(A8) 



+ k2 + ki logt + p2 log t 



dt 
T 



Finally, the scries (4.3) is obtained by subtracting the series (A6) from (A4). The resulting 
coefficients (4.4) arc calculated by means of (A5) and (A7): qi = /Ci+i — li+i for i — 
1, 2, . . . ,7. Since /Ci — Xi = 0, the contribution oc s~^logs is absent in (4.3). 
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Appendix B 

Now let us estimate /jr(p) from (4.6). First, we expand the products t Jo{N't)G{Kt) and 
tYo{J\ft)G{Kt) for a small argument t ^ 1: 

t JoiAft) G{Kt) ~ %t-^ +t (q2 log' {Aft) 

2 3 (Bl) 

+ 5^n3_, \og' [Mt)) + e J2^r-i log^Aft) + 



i=0 i=0 



tYo{Aft)G{Kt) ~ log(|AAt) + t (q2 V(AAt) 

3 4 (B2) 

i=0 i=0 

where the parameters ni,n2, ■ ■ ■ are expressed by means of the previously found coeffi- 
cients (4.4): 

ni = 3 log(|-^) g2 + 53, 

712 = 3 log^(;^) ?2 + 2 log (^-^) 53 + 54, 



Us = 



ris = - ^ + «^ + 3 log(^) ^4 ■ 

The coefficients mi,m2, ■ ■ ■ , ttiq are expressed by means of rzi, rz2, . . . , rzs and gi and 52: 
7711 = ^^§(2) ^2 + "-1 , m2 = ^°g(^) "-1 + "-2 , 

77l3 = log(l) ^2 + ri3 , 7714 = ^^+ log(|) ^3 , (^4) 

77i5 = 774 , 7776 = A/"^ ^ + ^°g(|) "-4 + "-5 ■ 

The coefficients 775 and 777, and 7777, 7773, and 7779 are present in (Bl) and (B2) only formally: 
their explicit values are not of importance for the present investigation. 

Using (B1)-(B4), we estimate the integrals that constitute It{p) and obtain 

I^p) ~ /o + /i log\Afp) + /2 log{Afp) 

3 (B5) 
+ p'J2 f-r-i log' i-^P) + P' h log' (AAp) + . . . , 

i=0 
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where 



f7^^ (Xy + (l - l0g|)x^) + log(l) 5^2-^3, 

The following notation is accepted in (B6): 



^2 

4 



(B6) 



„ _ n4 

/9 - -^g 



4 2 ' 2 2' 



In addition, we denote 



Ij = Ij- ^logAA, 



oo 1 

Tj = ^lim y uNt)G{Kt)tdt 

oc 1 

= J Yo{Aft)G{Kt)tdt + loggATt)^ 



where 



(B9) 



(BIO) 



It should be pointed out that possible contributions oc log^(A/'p) are cancelled in (B5): for 
instance, the coincidence of the coefficients 714 and ms (see (B4)) results in the absence of the 
corresponding term of the fourth degree in p. Moreover, (B3) and me (B4) are necessary 
to calculate fg (B6). 

Eventually, the estimate (4.7) appears after re-arrangement of (B5). Generally, the cor- 
responding coefficients Tq, ri, . . . , look as follows: 

ro = /o + fi log^^f + /2 logAA, ri = 2/1 logM + h , 

r2 ^ fi: ^ fi, r4 = 3 /4 log + /s , 

rs = 3/4 + 2/5 logAA + /e, (Bll) 

re = UWM + /slog'AA + /e logAA + A, 

r. = /. = (.^-AA^)^. 
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We obtain ro,ri,r2, and rs by means of (B6). To obtain r4 (4.8), we re-express f^: 

Now we turn to and rg from (5.10). Although /g and fj are left undone, their depen- 
dence on const can be extracted by means of (B6)-(B10) as follows: 

/6 = + /r = K'{q + logAA)^ + ... , (B13) 

where C and g are determined by (5.11). Since (B6) and /s (B12) are free from const, 
the equation (5.10) is indeed valid. 
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